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We calculate the high-temperature series of the magnetic susceptibility and the second and fourth 
QQ ' moments of the correlation function for the XY model on the square lattice to order f}"^^ by applying 

(^ I the improved algorithm of the finite lattice method. The long series allow us to estimate the 

("^ , inverse critical temperature as /3c = 1.1200(1), which is consistent with the most precise value given 

^Nj ' previously by the Monte Carlo simulation. The critical exponent for the multiplicative logarithmic 

correction is evaluated to be S = 0.054(10), which is consistent with the renormalization group 

prediction oi 6 = j^. 
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I. INTRODUCTION 



It is believed that the XY model in two dimensions exhibits a phase transition of the Kosterhtz-Thouless (K-T) 
type[l|, which is driven by the condensation of vortices. From the renormalization group arguments, it was predicted 



^ ■ that the correlation length has an essential singularity at the transition temperature T^ as 



?-exp(^-j, (1) 

(-H . with cr = ^ where t = T/T^ — 1 is the reduced temperature and 5 is a non- universal constant. At the critical 
Q ' temperature the correlation function behaves like 
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[~^ ' with f] — 1 ^-iid ^ ~ iki ''^^^^ the j'-th moment rrij of the correlation function behaves like 
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^^ ' The free energy and its temperature derivatives (i.e., the internal energy and the specific heat) were also predicted to 
OO , behave like 

o 
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f ^ ^ + non-singular term . (4) 



At the critical temperature, the first term on the right-hand side of Eq. ([¥]) has an essential singularity with itself and 
5^ [ its derivatives going to zero , while the second term stays nonzero. 

The behavior in Eq. ([T|) for the correlation length has been well established both by numerical simulations and the 
high-temperature expansion. The standard Monte Carlo simulation 0, Q gave f3c = 1.130(15) with 6 = 2.15(10) and 
Pc = 1.118(5) with b = 1.70(20) for the square lattice. (Here (3c is the critical inverse temperature, which will be 
defined below.) More precise values Pa = 1.1208(2) with b = 1.800(2) [| and Pc = 1.1199(1) with b = 1.776(4) g were 
obtained using the finite-size scaling technique and the renormalization group finite-size scaling method, respectively. 
In the latter approach, the renormalization group flow of the observable was matched with that of the exactly solvable 
BCSOS model. The latter value of Pc was recently confirmed by a large scale Monte Carlo simulation on a 2048 x 2048 
lattice using the finite-size scaling method [6|. On the other hand the high-temperature expansion for the magnetic 
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susceptibility and higher moments of the correlation function to order /3^^ [T4I and /3^^[TJ gave less precise value 
j3c = 1.118(3) with b = 1.67(4) and Pc = 1.1198(14) with b = 1.77(1), respectively. It seems that the available high- 
temperature series is not long enough to give the estimation of the values for the critical temperature and other critical 
parameters in the same precision as the Monte Carlo simulation. So it is desirable to extend the high-temperature 
series to much higher order. 

As for the critical exponent 9 for the multiplicative logarithmic correction in the moments of the correlation function, 
there has been controversial arguments. Negative values ranging from 9 = —0.077 to 9 = —0.056 0, H, Q were given 
by the analysis of the numerical simulation based on the thermal scaling formula ([3]), while the finite size scaling 
analysis gavepositive values in the range of 6* = 0.02 — 0.035@,[3, [l3,[li|. The large scale Monte Carlo simulation by 
Hasenbusch [<y\ gave a value 9 = 0.056(7), which is consistent with the renormalization group prediction, assuming a 
modified finite-size-scaling behavior 

mo -L2-"(C + In L)^^ (5) 

where L is the size of the lattice used in the simulation and C is a constant. On the other hand, the high-temperature 
expansion series to order /3^^ [i2| gave negative values of 6* = —0.042(5) and 9 = —0.05(2) assuming the thermal 
scaling ([3]). We should note that these values from the high-temperature series are obtained only for the Dlog-Pade 
approximants and the general inhomogeneous differential approximants do not give convergent result to this order. 
Higher order series would be needed again to resolve the discrepancy between the result of the high-temperature series 
analysis and the renormalization group prediction. 

A commonly known method for series expansions is the graphical method[15]. However, in this method, one must 
list all the graphs that contribute to the desired order of the series. An alternative and powerful method to generate the 
expansion series is the finite lattice method[la,ll3j[3. It avoids listing up all the graphs and it reduces the problem to 
the calculation of the partition functions for the relevant finite-size lattices, which is a rather straightforward procedure 
if we use the transfer matrix formulation. In many cases, the finite lattice method generates longer series than the 
graphical method [1^, [20, lU, 123, 123, H, [2^, l26|. Unfortunately, in the case of the XY model in two dimensions, the 
original finite lattice method can generate a high-temperature series that is at most as long as the series that can be 
obtained by the graphical method. 

In order to generate long series for the the XY model in two dimensions, we here apply an improved algorithm of 
the finite lattice method developed by the author and Tabata [23, l2a| ■ This improved algorithm is powerful in the 
case of models in which the spin variable at each site takes more than two values, including an infinite number of 
values. This algorithm was applied to generate a low-temperature series for the solid-on-solid model[27] and high- and 
low-temperature series for the g-state Potts model in two dimensions [23|. In both cases, it generates much longer series 
than the original finite lattice method. The XY model in two dimensions can be mapped to a kind of solid-on-solid 
model, and the improved algorithm of the finite lattice method in fact enabled us to obtain the high-temperature 
series for the free energy of this model on the square lattice to order /3^*[2^, which is two times longer than the 
series previously derived. From the analysis of the obtained long series we confirmed that the free energy of the two- 
dimensional XY model behaves like Eq. ([4]) , with values of the critical temperature and the non- universal constant 
b that are close to the values obtained in the study of the correlation length. We apply this improved algorithm of 
the finite lattice method to generate the high-temperature series to order (3^^ for the magnetic susceptibility and the 
second and fourth moments of the correlation function. The obtained long series will provide the value of the critical 
temperature in the same precision as the latest large scale Monte Carlo simulations and the value of the critical 
exponent 9 which is consistent with the renormalization group prediction. 

In section 2, we describe how to apply the improved algorithm of the finite lattice method to generate the high- 
temperature series for the moments of the correlation function. In section 3, the high-temperature series to order /S^^ 
are given. Section 4 is devoted to the analysis of the obtained series to evaluate the critical parameters. 

II. ALGORITHM 

We consider the XY model defined on the square lattice. The Hamiltonian of this system is 

-ff = - ^ JsiSj , (6) 

where Si is a two-dimensional unit vector located at the lattice site i, and the summation is taken over all pairs {i,j) 
of nearest neighbor sites. The correlation function is given by 

< Sr^so >= —— , (7) 



where Z is the partition function 



and 



r^.o = / ]^rf6'iSj;Soexp \-Tf]- 



(9) 



Here T is the temperature and di is the angle variable of the spin si — {cos 6i, sin^i). 

This model can be mapped exactly to a solid-on-solid model and the numerator and denominator of Eq. ([9]) can be 
rewritten as 

{h I — oo<- hi <.-\-oc} {iyj) 

and 

r.,o= E n^t''.-''.+A..;.,oi(/3)> (11) 

{/l I — OC</li< + Oo} {i:j) 

where (3 = J /kT, In is the modified Bessel function, the product is taken with respect to all the pairs of neighboring 
plaquettes, and the variable hi at each plaquette i takes integer value ranging from — cxd to +cxd. In Eq. pip the /S.ij;x.Q 
is defined as 

A _ J 1 for bi^j G L^fi, , , 

^'^r,x,o - I otherwise, ^ ' 

where 6^^ is the bond sandwiched by the neighboring plaquettes i and j, and L^.o is the set of bonds on an arbitrary 
shortest path connecting the sites — (0, 0) and x — (xi, X2) along the bonds, for which we adopt here the path that 
starts from the site and go straight first in 1 direction and then in 2 direction reaching the site x when a;ia;2 > 0, 
and first in 2 direction and then in 1 direction when xia;2 < 0. 

The improved algorithm of the finite lattice method to generate the high-temperature expansion series for the 
correlation function of this model is essentially the same as that for the free energy described in Ref. [23| ■ We first 
calculate the correlation function for each of the finite-size rectangular lattice A(/i x ^2,^) with a restricted range of 
the value of the plaquette variable 

Here the finite size lattice A(Zi xl2,p) is specified by its size |A| — li XI2 and its position p, and each plaquette variable 
hi is restricted so that h- < hi < h+ (where /i_ < and /i+ > 0) both in the calculation of the numerator and the 
denominator. We define the size of the finite lattice so that the li x I2 lattice involves li x I2 plaquettes, including the 
bonds and sites on their boundary. For instance, the 1x1 lattice consists of a single plaquette, including 4 bonds and 
4 sites. We take into account finite-size lattices with li = and/or I2 — 0. An ?i x lattice consists of li bonds and 
li + 1 sites with no plaquette. The 0x0 lattice consists only of one site with no bond or plaquette. The boundary 
condition is taken such that all the plaquette variables outside the li x I2 lattice are fixed to zero. The numerator and 
the denominator of Eq. (I13p can be calculated efficiently by the transfer matrix method using a procedure in which a 
finite-size lattice is built one plaquette at a time[30|, |3l|. 

We then define 0x,o(A, /i+, /i-) of the finite size lattice A and of the restricted range of the plaquette variables 
recursively as 

4>Xfi{^, h+, h-) =< SxSQ >A;/l+Jl_ 

0,,o(A',/i',,/i'_}. 



E 



A'CA, 0</i^</i_|_, h_<h'_<Q 

(14) 



It should be noted that 02;.o(A, /i+, /i_) = if the site a; or is not included in A. 
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FIG. 1: (a) Example of the main polymer in the standard cluster expansion. The open circle indicates the site (0,0) where 
the spin so exists and the closed circle indicates the site x — (a;i, 0:2) where the spin s^ exists, (b) Example of the sub-polymer. 
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FIG. 2: Examples of the cluster composed of a single main polymer that contributes to the lowest order term of the high- 
temperature expansion of (/>(A, h+,h-) for h+ > 1 and h- < —1. The size of the finite lattices in these examples is |A| = 5 x 5. 
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FIG. 3: Examples of the cluster composed of a single main polymer that contributes to the lowest order term of the high- 
temperature expansion of cf>{A,h+,h-) (a) for h+ > 1 and h- — and (b) for h+ — and h- < —1. The size of the finite 
lattice in these examples is |A| = 5 x 5. 



The correlation function in the thermodynamic limit is then given by 



< S^j:Sq > 



lim 

|A| — ^oo, /i+ — *oo, h^ 



< SxSo >A-h+,h. 



A,h+,/i_ 



(j)xfi{K,h+,h^). 



(15) 



In the last line of Eq. (fTS]) the summation should be taken for all the lattice size and all of its position and all integer 
values of /i+ and /i_ with < /i+ < 00, —00 < /i_ < 0. 

In the standard (graphical) cluster expansion of the correlation function for the SOS model, a cluster is composed 
of polymers: one main polymer that consists of the set of bonds Lxfi and connected plaquettes attaching to it and 
possible sub-polymers that consists of connected plaquettes. An example of the main polymer and sub-polymer can 
be seen in Fig. 1. A value hi [^ 0) is assigned to each plaquette i of the polymer. We can assign to each cluster two 
numbers, ft-max (> 0) and /imin (< 0), which are the maximum and the minimum, respectively, of the plaquette variable 
hi in all the plaquettes of the polymers of which the cluster consists. Then, we can prove[l7j that (j)x,o{A.,h+,h-) 
includes the contributions to < s^sq > from all the clusters of polymers in the standard cluster expansion for which 
^max = ^+ and /i,nin = /i- and that can be embedded into the lattice A but cannot be embedded into any of its 
rectangular sub-lattices. 



Now we consider from what order the series expansion of (j>x.o{-^, h+, /i_) with respect to P starts. It is enough to 
give the order for the position x — (xi, a;2) so that Xi > and X2 > 0. The order for the other cases is known by the 
90° rotational symmetry of the model. Any cluster that contributes to the lowest-order term of the series expansion 
for (j},j;fi{A, /i+, /i_) consists only of a main polymer. Hence the series expansion of (j)xfl{A, /i+, h^) begins from order 
pn^,o{A,h+,h-) in the case of /i+ > 1 and h^ < -1 with 



nx,oi^,h+,h_) = < 



2/1 + 2I2 - a;i - a;2 + 4/i+ + 4|/i_ | - 4 

for {pi,P2) = (0,0) and {p[,P2) = (a;i,a;2), 

2/1 + 2/2 - xi -X2+ 4/i+ + 4|/i_ I - 6 

for p2 = 0, p'l = xi, and (pijPa) ^ (0, X2) 

2li + 2/2 - xi ~ X2+ 4:h+ + 4|/i_ I - 8 
for the others. 



(16) 



Here we denote {pi,P2) as the position of the bottom-left corner and (^'1,^2) ^^ the position of the top-right corner, 
respectively, of the lattice A with its size h x I2 {p'l — pi = h and p'2 — P2 = h)- Examples of the main polymer are 
given in Fig. 2 (a) (b) (c), which correspond to the three cases in Eq. (fT6|) . respectively. In the case of /i+ > 1 and 
/i„ = we have 

2/1 + 2^2 - xi - 2:2 + 4/1+ 

for pi = 0, P2 = X2 and {p2,p'i) ^ (0,xi), 
2/1 + 2^2 - a;i - a;2 + 4/i+ - 2 

for pi = 0, P2 < and ( p^ > X2 or p[ = xi). 



na;,o(A, /i+,/i-) = < 



2/1 + 2I2 -x-Y- X2+ 4/1+ - 2 

for p'2 = X2, p'l > xi 

and ( pi < or (pi,p2) = (0,0)), 
2li + 2I2 - a;i - 2:2 + 4/i+ - 4 

for the others. 



(17) 



and in the case of /i+ = and /i_ < — 1 



nx,o{A,h+,h^) = < 



2?i -1-2/2 -xi -X2 +4|/i_ 

for {pi,p2) = (0,0) and (^,^2) = (a;i,a;2), 
2/1-1-2/2 -.Ti -X2-(-4|/i_|-h4 

for p2 = 0, p'l = xi and (pi,P2) ^ (0, 0:2), 

2/1-^2/2 -.Ti -X2-H4l/l_ 

for p2 = 0, pi < or p'^ = xi, P2 > 2^2, 

2/1 -1-2/2 -xi -X2 +4|/i_| -4 
for the others. 



(18) 



Examples of the main polymer are given in Fig. 3 (a) and (b), which correspond to the last case in Eg. p?)) and the 
last case in Eq. ([18]) , respectively. In the case of /i+ = and h^ — we have 



,. , x_ , /1+/2 for (pi,p2) = (0,0) and (p'i,p'2) = (a;i,X2), 



(19) 



Here n^fiiA, /i+, /i_) = 00 implies that (j)xfi{A, /i+, /i_) — 0. 

Thus, in order to obtain the expansion series for the correlation function < SxSq > to order f]-'^ , we have only to 
take into account all combinations of the rectangular lattice A and the range of the plaquette variable (/i+, /i-) that 
satisfy the relation nxfl{A., h+, /i_) < N in the summation of Eq. (fTS]) and to evaluate each of the (px,o{-A, h+, /i_) to 
order f3^ . 

The j-th moment of the correlation function is given by 



E 



Xy < SxSq > 



(20) 



where |a;| is the distance between the site x and 0. The moment can be calculated more efficiently in the following 
way. First we calculate Tj(li, /2, /i+, /i-) defined by 



rj(/i, /2, /i+, /i- 



^ \x-y\^rx^y{A,h+,h^) 

x.yQA 

— ZiA) 



(21) 



TABLE I: Coefficients of tfie high-temperature series for the magnetic susceptibility mo- 



n 




ai"^ 
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/ 
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1 
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1 


2 


12 


/ 


1 


3 


34 


/ 


1 


4 


88 


/ 


1 


5 


658 


/ 


3 


6 


529 


/ 


1 


7 


14933 


/ 


12 


8 


5737 


/ 


2 


9 


389393 


/ 


60 


10 


2608499 


/ 


180 


11 


3834323 


/ 


120 


12 


1254799 


/ 


18 


13 


84375807 


/ 


560 


14 


6511729891 


/ 


20160 


15 


66498259799 


/ 


96768 


16 


1054178743699 


/ 


725760 


17 


39863505993331 


/ 


13063680 


18 


19830277603399 


/ 


3110400 


19 


8656980509809027 


/ 


653184000 


20 


2985467351081077 


/ 


108864000 


21 


811927408684296587 


/ 


14370048000 


22 


399888050180302157 


/ 


3448811520 


23 


245277792666205990697 


/ 


1034643456000 


24 


83292382577873288741 


/ 


172440576000 


25 


376988970189597090587 


/ 


384296140800 


26 


62337378385915430773643 


/ 


31384184832000 


27 


480555032864478422139959 


/ 


119830523904000 


28 


80636088313579215330647 


/ 


9985876992000 


29 


19240186846097940775812460721 


/ 


1186322186649600000 


30 


34266867760374182809422566317 


/ 


1054508610355200000 


31 


9864232002328615891762221069959 


/ 


151849239891148800000 


32 


93697376428024822547085555709 


/ 


723091618529280000 


33 


665861878626519700317398249291449 


/ 


2581437078149529600000 



for each finite size lattice A. We note that the numerator and the denominator of Eg. lpTj) depend only on the lattice 
size |A| = li X I2 and on /14., ft,_ and that they are independent of the position of the lattice after the summation is 
taken for x and y. They can be calculated efficiently again by the transfer matrix method using the procedure in 
which a finite-size lattice is built one plaquette at a time not only for j = but also for j > 2 (j : even), the detail 
of which was described in the application of the finite lattice method to the calculation of the low-temperature series 
for the second moment of the correlation function in the simple cubic f sing model |23j . The moment is then given by 



Too 



E 



4'j{h,l2,h+,h^) 



(22) 



where 4>j{h, h, h+. ft,_) is defined recursively as 



E {ii-i[ + m2-i'2 + ^)hii'iA,K^h'_). 



(23) 






From Eq.dTH) - ^ we know that (j>j{li,l2, h+, h^) starts from order j3f^(h,i2:h+,h^) ^^^^ 

{h + I2 for h+ = h^ = 0, 
h+l2 + 4ft.+ - 4 for /i+ > and h^ = 0, 
h + I2 + ilh^l - 3 for /i+ = and /i_ < 0, 
I1+I2 + 4:h+ + 4\h_\ - 7 for ft.+ > and h_ < 0. 



(24) 



III. SERIES 



We have calculated the high-temperature expansion series to order /3 for the j'-th moments rrij (j = 0, 2, 4) of 
the correlation function for the XY model on the square lattice (to-q is the magnetic susceptibility). The obtained 



TABLE II; Coefficients of the high-temperature series for the second moment rn,2 of the correlation function. 



n 




a<=' 
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32 
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1 


3 


162 


/ 


1 


4 


672 


/ 


1 


5 


7378 


/ 


3 


6 


24772 


/ 


3 


7 


312149 


/ 


12 


8 


77996 


/ 


1 


9 


13484753 


/ 


60 


10 


28201211 


/ 


45 


11 


611969977 


/ 


360 


12 


202640986 


/ 


45 


13 


58900571047 


/ 


5040 


14 


3336209179 


/ 


112 


15 


1721567587879 


/ 


23040 


16 


16763079262169 


/ 


90720 


17 


5893118865913171 


/ 


13063680 


18 


17775777329026559 


/ 


16329600 


19 


1697692411053976387 


/ 


653184000 


20 


41816028466101527 


/ 


6804000 


21 


206973837048951639371 


/ 


14370048000 


22 


721617681295019782781 


/ 


21555072000 


23 


79897272060888843617033 


/ 


1034643456000 


24 


2287397511857949924319 


/ 


12933043200 


25 


5412508223507386985733313 


/ 


13450364928000 


26 


7139182711315236460182251 


/ 


7846046208000 


27 


107851995064346070336358789 


/ 


52725430517760 


28 


75355895214528595226953733 


/ 


16476697036800 


29 


1724076192091313972941261252343 


/ 


169474598092800000 


30 


53429382179619216000735913825117 


/ 


2372644373299200000 


31 


7534609247991680570113055733178247 


/ 


151849239891148800000 


32 


344373363823985360326961701207501 


/ 


3163525831065600000 


33 


68217033997582723452684940622596049 


/ 


286826342016614400000 



expansion coefficients are listed in Table I, II and III, where the coefficient a„ is defined as 



nii 



N 



,U) 



PY 



(25) 



We have checked that each of the (j)j{li,l2;h^,h^)'s in Eq. ([^5]) starts from the correct order in (3, as given by 
Eq. ((24)) . Our series coincide exactly with the series to order 0^^ for j = 0,2 and 4 given by Campostrini et al.|12| 
and to order (3'^^ for j = and 2 given by Butera and Pernici[13lll4|. which was obtained by a graphical method and 
a non-graphical recursive algorithm based on the Schwinger-Dyson equations, respectively. 

To obtain the large numerators and denominators of the coefficients exactly, we have used the same technique as 
in the calculation of the free energy series^]. In each step of the calculation the series expansion of a function /(/3) 
is expressed as 



/(«-i:S(f) . m-Ji^il 



then the product of the two functions is given by 

/(/3)ff(/3) = E 

n 

with 

n 

c« = E ~r 

-^ — ^ ri' 



n\ \ 2 



n'\{n — n')\ 



^n'^n — n' ? 



(26) 



(27) 



(28) 



and, if a„'s and 6„'s are integers, c„'s are also integers. 

The calculations were carried out on a PC cluster at the Information Processing Center at OPCT and on an 
Altix3700 BX2 at Yukawa Institute of Kyoto University. 



TABLE III: Coefficients of tfie high-temperature series for the fourth moment nn of the correlation function. 



n 




a<^' 










/ 




1 


4 


/ 




2 


96 


/ 




3 


930 


/ 




4 


6112 


/ 




5 


96850 


/ 




6 


147648 


/ 




7 


7305173 


/ 


12 


8 


2319540 


/ 




9 


498173873 


/ 


60 


10 


1271029508 


/ 


45 


11 


2210163319 


/ 


24 


12 


2606525954 


/ 


9 


13 


4449953438647 


/ 


5040 


14 


3300804041221 


/ 


1260 


15 


81597010130527 


/ 


10752 


16 


1952458419167723 


/ 


90720 


17 


782261299458533011 


/ 


13063680 


18 


222847854860540393 


/ 


1360800 


19 


287995465331747559427 


/ 


653184000 


20 


15925021359550756357 


/ 


13608000 


21 


8810941751514830879983 


/ 


2874009600 


22 


8551911058786623741001 


/ 


1077753600 


23 


21013548184238208070598633 


/ 


1034643456000 


24 


1108989303528609120610577 


/ 


21555072000 


25 


31566975982648926005193559 


/ 


244552089600 


26 


629546809875204592961945533 


/ 


1961511552000 


27 


1043070737244084798242371813021 


/ 


1318135762944000 


28 


10631253805062761391056806729 


/ 


5492232345600 


29 


5575566587216452758709818963741041 


/ 


1186322186649600000 


30 


2240485100208335200628893992637883 


/ 


197720364441600000 


31 


63403925202711962924615730240227527 


/ 


2336142152171520000 


32 


8174296006615548272296498779225887 


/ 


126541033242624000 


33 


394547895932847211470169206205154958361 


/ 


2581437078149529600000 



IV. SERIES ANALYSIS 



From Eq.Q the logarithm of the moment of the correlation function is expected to behave near the critical 
temperature like 



In) 



(2+J-r/)6 



+ 0(lnT), 



(29) 



where t — 1 — (3 / [3c and f3c — J /kTc is the critical inverse temperature. Here we analyze it by the first order 
inhomogencous differential approximation(IDA), in which the differential equation for /(/9) = \nnij is satisfied as 

g™(/3)/'(/3) + PmiiP) + RkiP) = 0(/3™+'+^-+^), (30) 

where QmiP), Pi{P) smd Rk{P) are polynomials of order to, I and k respectively and Qm(0) = 1. The critical inverse 
temperature Pc is given by the zero of Qm{P) and the exponent a is evaluated by 



PliPc) 



(31) 



In the analysis by the first order IDA here and below, we restrict tti + / + fc + 2 to be the maximum order of the 
analyzed series with —1 < fc < 9 and |?7i — Z| < 4, and exclude the approximants that have another zero of QmiP) 
with 1/3 — (3c\l I3c < 0.10, which is called near-by singularity. 

In Fig. 4 we plot the zeros of Q^ {(3) in the complex plane of /3. Besides clear accumulation of the points around 
f3 — 1.12 we see another accumulation around the nonphysical point (3 = (3o ^ —1.2 ± O.li. In order to remove the 
influence of this nonphysical critical point, we have made Euler transformation 



P' 



(3 



1 - (3/Po 



(32) 



with /3o = —1.2. In fact after this transformation the IDA gives a series of (/3c, <j) which converges better onto a 
straight line. Biased analysis fixing cr = ^ gives (3c = 1.1200(4). Hence we apply this Euler transformation in all of 
the series analysis presented below. 
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FIG. 4: Plot of poles in the unbiased inhomogeneous deferential approximation for In mo. The open circle indicates the poles 
corresponding to the physical singularity . 
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FIG. 5: /3c given by the biased IDA for In mo + cln 



4/3 ■ 



The above analysis ignores the existence of the sub-leading correction terms to the leading power-law singularity of 
In rrij in Eq. (j29p . So we have analyzed a combination of mo and m2 as 



InTTiQ 



, rn2 
c In , 

4/3' 



(33) 



searching the parameter c which will give more convergent result, anticipating that the sub- leading terms will be 
cancelled between the first and second terms in Ea. ((33)l . In Fig. 5 we plot the estimated value of Pc versus c in the 
biased analysis fixing a — ^. We find that the choice of the parameter c = 0.08 — 0.12 gives nicely converging result 
of (3c = 1.12007(4). We have also analyzed another combination as 



In mo + c' In (1 + c" 7712) 



(34) 



where c" is an arbitrary parameter which should be chosen so that the combination will give the most convergent result. 
We have obtained the most convergent result of /?c = 1.12003(5) for c' = 0.05 and c" = 0.58 and /3c = 1.11997(5) for 
c' = 0.03 and c" = 0.68. Thus we can safely estimate (3c = 1.1200(1). This value is quite consistent with the most 
precise value /3c = 1.1199(1) obtained by Hasenbusch from the large-scale Monte Carlo simulation ^6*. Our value is 
much more precise than the value (3c — 1.1198(14) by Butera and Pernici^lSi, J^] from the high-temperature series to 
order /3^^ . 

Assuming the critical behavior of Eq.([T]) and Q we can estimate the non-universal parameter b by the Fade 
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FIG. 6: Plot of the estimated value of 9 by the 1st order IDA for the test function versus B. 



approximation of 






In 1 



m2 
Too 



and 



A 

/3c 



In 1 



, TO4 
TO2 



b (35) 

6. (36) 

By searching the parameter c or c' that gives the most convergent estimation of b keeping Pc = 1.12000 and ct = — , 

both of the two give the same result of 5 = 1.758(1) for c — 3.15 — 3.23 and for c' = 0.806 — 0.812 respectively. Here 
we have used all of [ttt., I] Pade approximants with m > 14 and I > 14. This value is a bit smaller than b — 1.800(2) [j| 
and b = 1.776(4) [5[ obtained in the Monte Carlo simulations. 

Unfortunately the long scries do not improve the estimation of the exponent t] so much. For instance, the Pade 
approximation of the quantities 



In 1 + c 



TO2 

Too 



and 



In Too 



\n{l + cJ^ 



In TO,o 



(37) 



V 



(38) 



give the most convergent result of r; = 0.256(2) for c = 0.69 and ry — 0.227(2) for c' = 0.60, respectively. The Pade 
approximation of 



A 

Pc 



In 1 



,, TO2 



rjb 



(39) 



gives rjb — 0.429(5) for c" = 0.83 — 0.89 and, if we combine this with the above result b — 1.758(1), gives t] = 0.244(3) 
As for the exponent 9, Dlog-Pade analysis of the 20th order series of the quantities 



mo 



m^^ 



-2a0 



TO2 

Too 



2-r, 



TO2 
Too 



6 — ?7 



{1 + 0(r'" In r)} 



(40) 
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FIG. 7: Plot of the estimated value of 6 by the 2nd order IDA for the test function versus B. 
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FIG. 8: Value of 9 for each of k estimated by the 2nd order IDA of the combined quantity. 

gave estimations 9 — —0.042(5) and 9 — — 0.05(2) [l3|, which have the sign opposite to the renormahzation group 
prediction oi 9 = jq- The situation does not change even if our long series are used in the Dlog-Pade analysis of these 
quantities. The series of the two quantities to 33rd order gives 9 — —0.019(1) and 9 — —0.015(9), respectively. The 
IDA with A; > gives rather convergent values within the same k but quite scattered values for different k's. 

One possible reason why IDA for these quantities give scattered values for different k may be the existence of the 
sub- leading logarithmically singular term in Ea. (j40p . which comes from the correction factor in Eq.([3]). In fact the 
sub-leading logarithmic singularity can strongly disturb the correct evaluation of the leading power low exponent 2a9 



if this exponent is as small as 
(3^^ for a test function 



16 



We plot in Fig. 6 the estimated value of 9 by IDA for the expansion series to order 



' {l + T^iA + Blnr)} 



(41) 



with 9 



:f^, (7 = i, and A = plotted versus B. We find that the estimated value of 9 is quite sensitive to the 
amplitude B of the logarithmic term, and although each approximant with the same k gives rather convergent result 
for any fixed value of B, the approximants with different k give the estimation of 9 far from each other. 
Thus we have evaluated 9 using the combination of mo , TO2 and m4 as 



mo"(H-c'm2)" (l + c"m4)' 



(42) 



with 



^(2 - 770) + a'(4 - 770) + a" (6 - %) = 



(43) 



12 

and 

a + a +a" = 1 (44) 

and r/o — \, which is also considered to behave hke Ea. (|¥T|) in general. We can however anticipate that, by taking 
the combination of the three quantities, the subleading logarithmic term may be cancelled if we choose appropriate 
values for the parameters a, c' and c" . We have used the biased 2nd order IDA; 



T 



'^ ..(/?)/" (/3) + tQ™,(/3)/'(/3) + PK/3)/(/3) + ^fc(/3) = 0(/3™^+"^+'+'=+3) (45) 



with Qm2(0) = 1- The exponent —2a9 can be evaluated by the solution 7 of 



Pc' ^ Pc 



7(7 - 1) ^^^^^^ - ^'^2:11^ + p, ip^) ^ 0. (46) 



Here we adopt Pc = 1-1200. The 2nd order IDA can make more precise evaluation of 9 than the 1st order IDA for 
the function like Eq. (|4T|l . Fig. 7 is the plot of the estimation of 6 by the 2nd order IDA for the 33rd order series of 
the test function (|4ip with (^ — js, (^ — ^, and A = 0.3. Of course the 2nd order IDA also gives different values of 
6 for each of fc if B 7^ 0, but if B is small enough it can present precise estimation for the exponent of the leading 
singularity. By the analysis of the real combined quantity (|42|) we have found that the estimated values of 6 converge 
to 6* = 0.050(15) for all the range of— l<fc<9ina domain of the set of parameters —a = 3.6 — 3.0, c' — 2.8 — 4.0 
and c" — 2.8 — 4.0. In the analysis we have restricted TO2 + mi + / + A; + 3 = 32 with —l<k<9 and 1112 > 7, mi > 7, 
I > 7. The approximants that have near-by singularity (i.e. the zero of Q^iP) or Qi{P)) with \P — Pd/Pc < 0.2 have 
been excluded. The most convergent result 9 = 0.054(10) is obtained for a = —3.143, c' = 3.134 and c" = 3.139. The 
values of 9 for each k in this set of parameters are shown in Fig. 8. We note that this value 9 — 0.054(10) is consistent 
with 9 = jt: predicted by the renormalization group. 

V. SUMMARY 

We calculated the high-temperature series for the zeroth moment (magnetic susceptibility) and the second and 
fourth moments of the correlation function in the XY model on the square lattice to order P^^ by using the improved 
algorithm of the finite lattice method. The obtained long series have presented us an estimation for the value of the 
critical inverse temperature as Pc = 1.1200(1), which is consistent with the most precise value given previously by 
the Monte Carlo simulation. The critical exponent 9 for the multiplicative logarithmic correction is evaluated using 
the combination of the three moments of the correlation function, giving 9 — 0.054(10), which is consistent with the 
value 9 = jq predicted by the renormalization group argument. 
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